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Abstract: We show the existence of a duality between the c-map space describing the 
universal hypermultiplet at tree level and the matrix model description of two-dimensional 
string theory compactified at a self-dual radius and perturbed by a sine-Liouville potential. 
It appears as a particular case of a general relation between the twistor description of four- 
dimensional quaternionic geometries and the Lax formalism for Toda hierarchy. Furthermore, 
we give an evidence that the instanton corrections to the c-map metric coming from NS5- 
branes can be encoded into the Baker- Akhiezer function of the integrable hierarchy. 
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1. Introduction 

Integrability is a very powerful tool to find exact solutions and it is always a great achievement 
when a theory is proven to possess an integrable structure. A particularly interesting class 
of theories where various appearances of integrability have been observed during the last few 
years is characterized by the presence of N = 2 supersymmetry. Remarkably, it appears both 
in the gauge |], 0, |3|, £|, [5| and in the string theory context || [7). However, an exhaustive 
explanation of these observations is missing so far. 

In this note we reveal one more instance of an integrable structure in theories with N = 2 
supersymmetry. Namely, we observe that the four- dimensional quaternion-Kahler (QK) space 
obtained by the local c-map construction || || is in a certain sense dual to a sine-Liouville 
perturbation of the c = 1 string theory compactified at the self-dual radius. The former 
represents the tree level moduli space of the so called universal hypermultiplet |TD[ arising 



in compactifications of type II strings on a rigid Calabi-Yau. On the other hand, the latter 
is known to be an exactly solvable theory which is described by Toda integrable hierarchy 



11, 12|. In particular, the lowest equation of this hierarchy is given by the Toda equation, 



which is also known to parametrize four- dimensional quaternionic metrics with one isometry 
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n|. The duality which we report here is based on the observation that the two systems are 
described by the same solution of the Toda equation. 

However, the duality goes much further as it allows to identify also various quantities 
characterizing either the c-map or the c = 1 string. These identifications originate in a 
general correspondence between a twistor description of quaternionic manifolds and the Lax 
formalism for integrable hierarchies [ 14 1 . The duality between the c-map and the c = 1 string 
is just its simplest example. 

The above mentioned correspondence is not at all surprising. First of all, the existence of 
integrable hierarchies governing the four-dimensional hyperkahler (HK) geometry is very well 
known (see, for instance, [|T5| , |TE| , |T7 ). Moreover, recently it was realized that hyperkahler 



manifolds with a rotational isometry endowed with a hyperholomorphic line bundle can be 
associated with QK manifolds of the same dimension also possessing an isometry [Lj| . 



In particular, this QK/HK correspondence identifies the rigid and local c-map spaces and 
makes a contact between QK geometry and integrability. Although, this might seem to 
imply that the duality we propose is rather trivial from the mathematical point of view, its 
physical realization is quite fascinating as it suggests that the hypermultiplet sector of N = 2 
superstring compactifications and non-critical c = 1 string theory are closely related to each 
other. We are not aware of any examples where some perturbation of two-dimensional string 
theory was related to the hypermultiplet moduli space. 

Of course, this duality would become really interesting if it can be extended beyond 
the classical limit and used to extract quantum corrections on one side from the known 
results on the other side. We leave the systematic study of this problem for the future. 
Nevertheless, in this note we make the first step in this direction. Namely, we observe that 
the holomorphic functions on the twistor space of the c-map, which have been shown to 
generate NS5-brane instanton corrections to the tree level universal hypermultiplet, can be 
interpreted in terms of the Baker- Akhiezer wave function of the Toda hierarchy. Together 
with previous findings || |7|], this suggests that the whole non-perturbative picture describing 
N = 2 string compactifications possesses an integrable structure where the relations discussed 
here will find their natural place. 

2. Toda lattice hierarchy 

We start from a brief review of two-dimensional Toda lattice hierarchy |2U| which provides the 
general framework for what we are going to discuss. This integrable hierarchy describes many 
physical phenomena which, in particular, include a certain class of deformations of c = 1 string 
theory. More precisely, it encompasses all perturbations generated by the spectrum of tachyon 



operators corresponding to the c = 1 string compactified on a circle |TT| [12]. Furthermore, the 
Lax formalism of this hierarchy turns out be closely related to the matrix model description 
of the c = 1 string, which is provided by Matrix Quantum Mechanics (MQM). Therefore, first 
we present this formalism, then we explain its interpretation from the MQM point of view, 
and finally we give a particular solution of the Toda hierarchy which describes c = 1 string 
theory with the so-called sine-Liouville perturbation. Although we do not need the Toda 
hierarchy in its full generality to explain the relation with the c-map spaces, it is nevertheless 
useful to recall it as we want to discuss further extensions and generalizations of the proposed 
duality. 
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2.1 Lax formalism 

Let us first give a formal definition of the Toda lattice hierarchy. To this end, introduce two 
semi-infinite series 

oo oo 

L = r(s)CC + ^u k (s)ou~ k , Z = w _1 r(s) + ^w fc w fe (s), (2.1) 

k=0 k=0 

where s is a discrete variable labeling the nodes of an infinite lattice, Cj = e h9 / 9s is the shift 
operator in s, and the Planck constant h plays the role of a spacing parameter. The operators 
(271) are called Lax operators. Their coefficients r, u k and u k are taken to be functions of 



two infinite sets of "times" {t±fc}^=i- To each time variable one associates a Hamiltonian H± k 
generating an evolution along t± k by the usual rule 



h WT k = i H -k, L ], h WT k = [H-k,L]. 



(2.2) 



This system represents the Toda hierarchy if the Hamiltonians satisfy an additional require- 
ment, namely that they can be expressed through the Lax operators ( |2.1| ) as follows 



H k = (L% + i(L fc ) , H. k = (L fc )< + l(L fc ) , (2.3) 

where the symbol ( ) > means the positive (negative) part of the series in the shift operator Cj 
and ( )o denotes the constant part. The integrability follows from the existence of an infinite 
set of commuting flows generated by h-^- — H k whose commutativity is a consequence of (|2.2|). 

This definition shows that the Toda hierarchy is a collection of non-linear equations of the 
finite-difference type in s and differential with respect to t k for the coefficients r(s,t), u k (s,t) 
and u k (s, t). Its hierarchic structure is reflected in the fact that one obtains a closed equation 
on the first coefficient r(s,t) and its solution provides the input information for the following 
equations. It is easy to show that r(s, t) should satisfy a discrete version of the Toda equation 

Ot\ut—\ 

which gives the name to the hierarchy. 



For our purposes it is important to introduce the following Orlov-Shulman operators [21 



M = ^2kt k L k + s + ^2v k L- k , M = — ^2 kt^ k L k + s — ^ v_ k L~ k . (2.5) 

k>l k>l k>\ k>\ 

The coefficients v± k (s, t) are supposed to be found from the condition on the commutators of 
M and M with the Lax operators 

[L, M] = hL, [L, M] = -hi. (2.6) 

The importance of the Orlov-Shulman operators comes from the fact that they can be con- 
sidered as perturbations of the simple operators of multiplication by the discrete variable s, 
whereas the Lax operators appear to be a "dressed" version of the shift operator Cj. Indeed, 
if one requires that v± k vanish when all t± k = 0, then in this limit M = M = s, L = Cj and 
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L = Cj 1 . Moreover, the commutation relations (p.6|) can be recognized as a dressed version 
of the trivial relation 

[u,s] = hu. (2.7) 

Another application of the Orlov-Shulman operators is that they lead directly to the 
notion of r-function. Indeed, one can show that their coefficients must satisfy |22| 

au-W k - (2 - 8) 

This implies that there exists a generating function r s [t] of v±j t such that 

Vk(s,t) = n 7- . (2.9) 

Otk 

It is called the r-function of Toda hierarchy. It encodes all information about a particular 
solution. For example, it allows to extract the first coefficient in the expansion of the Lax 
operators 

Plugging this relation into ( p.4|) , one obtains that the r-function is subject to the Toda 
equation. In physical applications the r-function usually coincides with partition functions 
and the coefficients v & are the one-point correlators of the operators generating the commuting 
flows iffc. 

Finally, note that the Toda hierarchy can be equivalently represented as an eigenvalue 
problem given by the following equations 

xV = LV, hx— = M^, h— = H k m, (2.11) 
ox ot k 

plus similar equations for \I/ where L, M should be replaced by L, —M. The eigenfunction 
\l/(x; s, t) is known as the Baker- Akhiezer function. Often it has a direct physical interpretation 
and plays a very important role. 

2.2 String equations 

The Toda hierarchy provides a universal description for many integrable systems because the 
equations of the hierarchy have many different solutions. There are two ways to select a 
particular solution. The most obvious one is to provide an initial condition which can be, for 
instance, the r-function for vanishing times corresponding physically to the partition function 
of a non-perturbed system. However, the Toda equations involve partial differential equations 
of high orders and require to know not only the r-function at t± k = but also its derivatives. 
Therefore, it is not always clear whether the non-perturbed function is sufficient to restore 
the full r-function. 

Fortunately, there is another way to select a unique solution of the Toda hierarchy. It 
relies on the use of the so-called string equations which are some algebraic relations between 
the Lax and Orlov-Shulman operators. The advantage of the string equations is that they 
represent, in a sense, already a partially integrated version of the hierarchy equations. For 
example, instead of solving a second order differential equation on the r-function, they reduce 
the problem to certain algebraic and first order differential equations. 
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It is important, however, that the string equations cannot be chosen arbitrary because 
they should be consistent with the commutation relations (|2.6|) . For example, if they are 
given by two equations of the following type 

L — f(L, M), M — g(L, M), (2.12) 

the operators defined by the functions / and g must satisfy 

[f(u,3),g(u,8)]=-hf. (2.13) 

2.3 Dispersionless limit 

In this paper we are mostly interested only in the classical limit of the Toda hierarchy, which 
is obtained by taking the spacing parameter h to zero. In this limit one recovers what is 
known as dispersionless Toda hierarchy |22[. It is described by the same equations as above 
with the only difference that all commutators should be replaced by Poisson brackets and all 
operators are considered now as functions on the phase space parametrized by s and u with 
the symplectic structure induced from ( |2.7| ) 

{lo,s} = u. (2.14) 

As in the full theory, a solution of the dispersionless Toda hierarchy is completely char- 
acterized by a dispersionless r-function. In fact, it is better to talk about free energy since it 
is the logarithm of the full r-function that can be represented as a series in h 

\ogT = Y,K~ 2+2n Fn. (2.15) 

n>0 

The dispersionless limit is described by the first coefficient of this series, the dispersionless 
free energy Fq. It satisfies the classical limit of the Toda equation 

dtA-rFo + e^^O, (2.16) 



and can be selected from all solutions by the same string equations (|2.12| ) as in the quantum 
case. 

The quasiclassical asymptotics of the Baker- Akhiezer function is given by the usual WKB 
approximation. Since \1/ diagonalizes the Lax operator, it can be considered as a wave function 
representing the quantum state of the system in the L-representation. The conjugate variable 
is provided by M which implies that 



1 



log X 



^~exp<|-/ MdhgL}, (2.17) 



where M is considered as a function of L, s and tk- Using ( |2.5| ), one finds 

= e * s+o(»°) S = Y t k x k + a logx + I - Y ^ x~ n . (2.18) 



n 

k>l k>l 



The integration constant <fi can be related to the free energy as <f){s,t) = —d s Fo- 
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2.4 Example: MQM with sine-Liouville perturbation 

Two-dimensional or non-critical c = 1 string theory is known to be an example of integrable 
model described by the Toda hierarchy. However, this is difficult to see in its CFT formulation. 
On the other hand, the integrability becomes transparent in its matrix model description in 
terms of Matrix Quantum Mechanics (see P3L B3 for a review). 



To see how it appears, one should restrict to the singlet sector of MQM where all angular 
degrees of freedom are integrated out. The remaining matrix eigenvalues describe a system of 
free one-dimensional fermions in the inverse oscillator potential. The target space description 
of string theory arises as an effective theory of collective excitations of the free fermions. In 
particular, different backgrounds correspond to different states of the Fermi sea: the simplest 
linear dilaton background is dual to the ground state, whereas various perturbed backgrounds 
with a non-trivial condensate of the tachyon field come from ceratin time-dependent states. 

The most convenient way to describe the dynamics of these free fermions is to use the 
light-cone coordinates \T2\ x ± = where (x,p) are the coordinates on the phase space of 



one fermion. It is clear that x ± satisfy the canonical commutation relations 

{x + ,x_} = l. (2.19) 

In the quasiclassical approximation the problem reduces to finding the exact form of the 
profile of the Fermi sea in the phase space, given its asymptotic form at x ± — > oo. Here we 
are interested in deformations generated by the spectrum of the theory compactified on a 
circle of radius R. Then the profile of the Fermi sea is determined by the compatibility of the 
following two equations 

x+ x_ = m ± ( X± ) = ±^Y1 kt ± k x i /R +f i± ^Yl v ± k x ± k/R - ( 2 - 2 °) 

k>l k>\ 

Here x k J R correspond to the tachyon vertex operators with momentum k/R, t±k are their 
coupling constants, /i is the Fermi level, and v±k are to be found as functions of t±k and 
fi. For vanishing coupling constants, the equations ( |2.2Up describe the hyperbola x + x_ = fi 



which is the trajectory of a free fermion in the inverse oscillator potential. 

Comparing the r.h.s. of ( |2.20| ) with (|2.5| ), one immediately sees that M± are related to 



the Orlov-Shulman operators. In fact, one can establish a precise correspondence between all 
data of the Lax formalism and MQM quantities: 

• the lattice parameter s is identified with the Fermi level //; 

• the Lax operators coincide with the light-cone coordinates, L = x^ R , L = x]/ R ; 

• the Orlov-Shulman operators encode the asymptotics of the Fermi sea and are given by 
M = R.M.. M = RM : 

• the Baker- Akhiezer function coincides with the perturbed one-fermion wave function; 

• the string equations are the same as the equations for the profile of the Fermi sea and 
correspond to the choice / = (M / R) l l R L~ l , g = M. 
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The commutation relation ( |2.19| ) then follows from a combination of ( |2.6| ) with the string 



equations. Finally, the classical limit of the shift operator u turns out to be a uniformization 
parameter for the Riemann surface described by ( |2 . 2 Q|) . 1 

Note that the two Baker- Akhiezer functions, $ and give one-fermion wave functions in 
the two chiral representations, x + and x_ , respectively. The requirement that they represent 
the same quantum state is precisely what generates the profile equations ( |2.20|) . Due to this 
condition, the two wave functions are related by a Fourier transform. 

Restricting to a finite set of non- vanishing couplings t±k, it is straightforward to find an 
exact solution, i.e. an explicit form of the Lax and Orlov-Shulman operators as well as the 
r-function of this system. We are interested in the particular case R = 1 and t±k = for 
all k > 2. The perturbation induced by the couplings t±i is known to correspond to the 
sine-Liouville operator on the string worldsheet, and the compactification radius R — 1 is 
the self-dual point with respect to T-duality which leads to an enhanced symmetry. The 
corresponding solution of the Toda hierarchy is given by Hl2 



F ° = \^ ( log/i ~ ^) "A** 1 *- 1 ' x ± = y/Jlu^ 1 Tt T i. (2.21) 
3. Four-dimensional QK spaces and integrability 

Now we turn to a seemingly completely different subject — quaternionic geometry, which is 
characterized by the existence of three (almost) complex structures J 1 satisfying the algebra 
of quaternions and appears in two forms, hyperkahler and quaternion-Kahler. In fact, it is 
well known that in four- dimensions HK spaces possess an integrable structure, which in the 
presence of an isometry is given by the dispersionless Toda hierarchy presented in the previous 
section. Moreover, as we explained in the introduction, the QK/HK correspondence allows 
to transfer these results to the realm of four-dimensional QK geometry. In this section we 
establish precise relations between the quantities used to describe QK spaces with one Killing 
vector and the ones entering the Lax formalism. 

3.1 Toda ansatz 

The most direct way to see the appearance of Toda in the description of QK spaces is to 
recall that in four dimensions QK manifolds coincide with self-dual Einstein spaces with a 
non- vanishing cosmological constant. Such spaces are known to be classified by solutions of a 
single non-linear differential equation | 25fl . In the presence of an isometry, it can be reduced 



to the three-dimensional continuous Toda equation 



d z d- z T + d 2 p e T = 0, (3.1) 



which is identical to the equation ([2.1 6| ) on the dispersionless free energy of Toda hierarchy. 



In terms of the solution T(p, z, z), the metric is given by |fj6| 

d4 = — [j 2 (dp 2 + Ae T dzdz) +j^ 2 ( de + ©) 



(3.2) 



1 In addition, to describe the theory in Lorentzian signature, the Planck constant should be continued 
analytically to imaginary values, fi/roda = ifiMQM- 
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where the isometry acts as a shift in the coordinate 6. Here P = 1 — ~ pd p T, A is the 
cosmo logical constant, and G is a one-form such that 

d6 = i(d z Pdz - dzPdz) A dp — 2i d p (Pe T )dz A dz. (3.3) 



The integrability condition for (|3.3|) follows from G3.1|) . Of course, not only the metric, but 
also all geometric data such as quaternionic two forms, Levi-Civita connection, etc., can be 
expressed through T. We refer to [[27] for their explicit expressions. 

Comparing the Toda equations ( |3.1| ) and ( |2.16| ), one immediately finds the following 
identifications 

z, z t±i, p s, To <9 2 P . (3.4) 

In principle, one can go further and find similar relations for other quantities. However, 
as we already know, working with differential equations is not the best way to describe a 
solution of the hierarchy. It is much better to use string equations which turn out to have a 
certain twistorial interpretation |15, 14] . Similarly, QK manifolds have a much more powerful 



description in terms of their twistor spaces |28| . As we will see now, the two become essentially 
identical. 

3.2 Twistor space and Lax formalism 

The twistor space Z of a QK manifold A4 is a CP 1 bundle over A4, whose connection is given 
by the SU(2) part p of the Levi-Civita connection 2 on Ad and the CP 1 fiber corresponds to 
the sphere of almost complex structures J\ The twistor space carries a so called complex 
contact structure |29|, [30|, which is represented by a holomorphic one- form X such that the 
associated top-form X A (dX) n is nowhere vanishing. Locally it can be written as 

X W = - e* H (dt + p + - ip 3 1 + p- t 2 ) , (3.5) 

where the index [i] labels open patches of a covering of Z, t parametrizes the CP 1 fiber, 
and the function ^(rc^t) is known as the contact potential. In the case we are interested 
in, when A4 features one Killing vector, the contact potential is actually real, globally well 
defined, and t-independent so that one can put $W = 0(x M ). 

Locally, it is always possible to choose such coordinates that the contact form X acquires 
a standard form. Restricting to the case dim^ A4 = 4, this means that in a patch Ui C A4 
one can write 

AfM =daM + £ [i] dfM, (3.6) 

where (£, £, a) are analogous to holomorphic Darboux coordinates of complex symplectic 
geometry. Generically, these coordinates must be regular in Ui. The only exception is the 
coordinates defined around the north ({t = 0} G U + ) and the south ({t = oo} G U-) poles of 
CP 1 . In these patches the Darboux coordinates are assumed to have the following expansions 



oo 



ew(t) = nt^ +]re ] t ± ", l [±] (t) = E^ ] t ±r 



' " " " (3 7) 



'OO 

a W(t) = ^aWt ± "-2clogt, 

ra=0 



2 Recall that the holonomy group of a 4n-dimensional QK manifold is contained in Sp(n) x SU(2). 
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where 1Z is a real function on Ai and c is a real number known as anomalous dimension. Thus, 
£M have a simple pole and have a logarithmic singularity controlled by c. The latter has 
an important physical meaning generating the one-loop correction to the local c-map 



On the overlap of two patches Ui DUj, the two coordinate systems are related by a com- 
plex contact transformation. Such transformations are generated by holomorphic functions 



which we call transition functions. Together with the anomalous dimen- 



sion, they encode all geometric information about Ai and its twistor space Z. To extract the 
metric from them, the most non-trivial step is to find the Darboux coordinates as functions 
of t and coordinates on Ai. Although this cannot be accomplished generically, it is possible 
to write certain integral equations on the Darboux coordinates in terms of H™. We refrain 

l2| for more details. 



from writing them explicitly and refer to |31 



This twistor description has been related to the one based on the Toda ansatz ( |3.2|) in 
27j . For our purposes it is sufficient to provide the relation between the coordinates, contact 



and Toda potentials, which read as follows 



P 



I tt+] 
2 ^o > 



6 = — Imal 



T = 2 log ft. 



(3i 



Note that the only effect of the anomalous dimension is to shift the contact potential intro- 
duced in (B.5I), i.e. = e 



■- e^| c= o — c. The identification ( |3.8| ) nicely maps this shift to the 
ambiguity in the solution of the Toda equation, T(p, z, z) = T(p + c, z, z). Although the two 
solutions, T and T, are related by a coordinate change, due to the explicit dependence of the 
metric (|3.2|) on p, the geometry of the corresponding QK manifold is strongly affected by the 
parameter c. 

After this very brief review of the twistor description of QK spaces, one can easily see 
how it is embedded into the Lax formalism of the dispersionless Toda hierarchy. Indeed, 
comparing the expansions ( |3.7|) with ( |2.1|) , one observes a close similarity between the Lax 
operators and the Darboux coordinates £W; they are both given by Laurent series which 
contain one singular term. This similarity is further strengthen by the comparison of the 
coefficients of these singular terms. Using fl3.8| ) and ( j2.10 ), one finds 

K = e T ' 2 e 9 ° F °/ 2 = r, (3.9) 

in agreement with ( p.4|) . Correspondingly, the phase space variable u should be identified 
with the CP 1 coordinate, u^ 1 -H- t. Besides, the symplectic structures derived from (|2.6|) 
and dX Q3.6|) suggest that the Orlov-Shulman operators should be related to the Darboux 
coordinates Finally, it is well known (see, for instance, []15[ p2[ ) that the string equations 
( |2.12| ) can be considered as gluing conditions between the two patches around the north and 
south poles of CP 1 . The condition ( |2.13| ) in turn requires that the gluing conditions generate a 
symplectomorphism preserving the symplectic structure. We summarize all the identifications 
in the following table: 



lattice variable s 
d 2 s F 

t±i 
L,L 
M/L, M/L 
string equations 



t 

p = e 4 > 
T = 2\ogK 

z, z 

f[±] 

gluing conditions 



(3.10) 
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These identifications are however a bit incomplete because we need also to match the 
subleading terms in the expansions of the Lax operators and Darboux coordinates. This 
depends on the precise matching of string equations with the gluing conditions and should 
be considered case by case. In particular, in the next section we consider the situation of our 
special interest where all identifications can be made very precise and explicit. 

In addition, they say nothing about the role on the Toda hierarchy side of the third Dar- 
boux coordinate a. The point is that, in the case with an isometry, this is a derived quantity, 
namely, it can be given in terms of the transition functions and the explicit expressions for 
£ and £. Therefore, it is not essential for defining the system. Nevertheless, we will see in 
section || that it also has a counterpart in the Lax formalism being related to the quasiclassical 
asymptotics of the Baker- Akhiezer function. 



4. c-map vs c = 1 string 

4.1 Universal hypermultiplet at tree level 

Let us now consider a very special QK space obtained by the local c-map construction H, |?| , 
which maps a special Kahler manifold characterized by a homogeneous holomorphic prepo- 
tential F(X) to a QK manifold. Physically, the local c-map produces the tree level metric 
on the hypermultiplet moduli space which describes the low energy effective action of type 
II A string theory compactified on a Calabi-Yau manifold 2). Here we are interested in the 
case where the corresponding quaternionic manifold is four-dimensional. In this case the lo- 
cal c-map gives the moduli space of one hypermultiplet known under the name of universal 
hypermultiplet, which appears in the case where 2) is a rigid Calabi-Yau, i.e. with vanishing 
Hodge number ft, 2 ,i(2)). 

Since in our case the prepotential, we start with, is a holomorphic function of one variable 
and homogeneous of degree 2, it is always quadratic, F(X) = f X 2 . Although the modulus r 
carries an important physical information [33], for simplicity we restrict to the case r = — |. 



Then the resulting quaternionic metric is given by (|3.2|) where the potential T should be 
chosen to be one of the most trivial solutions of the Toda equation [R4L K5[ 



T = logp. (4.1) 

The coordinates (p — e^, z, 6) describe respectively the dilaton, a pair of RR-fields, and the 
NS-axion dual in four dimensions to the antisymmetric B field. 

On the other hand, the twistor description of this space is based on the covering by two 
patches only: U + covers the north hemisphere of CP 1 and U- covers the south one. The 
single transition function is determined by the prepotential, 3 

H [+ ~ ] = -4 Im = £ 2 . (4.2) 

In this case the Darboux coordinates on the twistor space have been found in [25] and, being 
expressed in terms of coordinates appearing in the Toda ansatz (|3.2|), read as follows 

£ [±1 = -(^ + ^) + v / p(t- 1 -t), 
£~W =2{z + yfpt), £ H = 2 (-z + ypt- 1 ) , (4.3) 
a W = 2i0 =f (p - 2( Re zf T 4^p( Re z)t ±x - pt ±2 ) . 
3 We adapted the normalizations to have simple duality relations. 
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4.2 Duality 

After all these preparations, it is now trivial to establish the correspondence between the 
four- dimensional c-map and the perturbed c = 1 string at the self-dual radius whose solution 
is encoded in ( |2.21| ). Indeed, since the Fermi level \i is identified with the coordinate p, both 
systems are described by the same solution of the Toda equation 

0jF o = T = logp. (4.4) 

Moreover, it is trivial to check that the identifications (|3.10|) imply 

x + =L= e +] + \ fW x_ = L = -£H + 1 £H 

! 2 x 2 (4-5) 

M/L =-f [+] , M/L = -t}-\ 



As it should be, the string equations (|2.20|) appear then as a consequence of the gluing 
condition 

|[+] _|[-] = _ 2 ^[±]. (4. 6 ) 

These equations provide the completion of the identifications in the general dictionary (|3.10|) 
and establish the full equivalence of the two systems. 

In fact, the relations (O) are not very illuminating. One can obtain much nicer identifi- 
cations if one introduces an additional patch on the twistor space of the c-map. Let it be an 
intermediate patch Uq around the equator of CP 1 (see Fig. |1|) and the associated transition 
functions be [ES| 

H [+o] = 2LF(£), fl-l-o] = 2LF(f). (4.7) 

This patch becomes especially important when one incorporates instanton corrections to the 
hypermultiplet moduli space [36|] and has an advantage that the Darboux coordinates in Uq 
transform nicely under symplectic rotations |J7]]. They are given explicitly by 

£ = p ] = (z - z) + y/p (t- 1 + t) , ( 4. 8) 
a = 2i ha® + e [01 | [0] ) = a + Ai^/p ^t" 1 + zt) , 

where we denoted o = —88 — 2i(z 2 — z 2 ). Applying ( |4.5|) , one finds remarkably simple relations 
between the Darboux coordinates (£,£) and the coordinates on the MQM phase space (x,p) 

i = x + — x_ = \/2p, £ = x + + x_ = \plx. (4.9) 



5. NS5-branes and the Baker- Akhiezer function 

In the previous section we clarified the meaning of the Darboux coordinates £ and £ under the 
correspondence with the c = 1 string. Does the third coordinate a have some interpretation 
on the dual side? To answer this question, let us combine ( p.6|) and (|3.5|) and integrate the 
resulting relation over a contour on CP 1 with a free upper limit. This leads to 

a = 2plogt + / £d£ + const, (5.1) 
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Figure 1: Covering of CP 1 by three patches, the associated transition functions, and the two 
contours incorporating the effects of NS5-branes. 



where the last term is independent of t. Then the dictionary ( |3.10| ) implies that the second 
term should be related to the phase of the quasiclassical Baker- Akhiezer function ( J2.17 ) . 

This expectation indeed turns out be correct and can be made quite precise in our particu- 
lar case characterized by the relations (|4.9| ). Since the string equations imply that M/L = x_ 
and M/L 



x + , the phases of the two Baker- Akhiezer functions, \I> and are given by 



S(x. 



x dx 



+ > 



S{x_ 



x^dx 



(5.2) 



Then using ( 2.21 ) and ( |4.8|) . it is easy to verify that 



-ia + £ 2 -£ 2 -plogt, 



o + S(x, 



--a + S(x_) = -- ia + r-C -plogt 



(5.3) 



This confirms the identification of a (or a combination thereof with other coordinates) with 
the WKB phase and shows that the coordinate a on the hypermultiplet moduli space can be 
seen as a common constant imaginary phase in the Baker- Akhiezer functions. 

Furthermore, the result ( |5.3| ) has a remarkable consequence. It is based on the observation 
that the terms in the brackets on the r.h.s. appear in the exponential of the holomorphic 
functions 

(5.4) 



which have been shown to incorporate NS5-brane instanton corrections to the classical ge- 
ometry of the universal hypermultiplet |]27|. Indeed, let C± be the two contours connecting 
t = (t = oo) to the point t±, corresponding to the complex submanifold £ =F £ = of the 
twistor space and found to be 



t. 



z 



t. 



y/P 

z 



(5.5) 



The constant term can be evaluated using (|2.1i 
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We depicted these contours in Fig. [I]. Then, being integrated over C±, the holomorphic 
functions -£^| 5 produce exponential corrections of the following form 

f — #NS5 ^ e -A7Tk(p+zz) T 7Tika _ ^ g-j 

This is a typical behavior of NS5-brane instantons since in our conventions p ~ g~ 2 and 



a is the axionic coupling associated with NS5-branes [[38]. The parameter /c is an integer 
interpreted as NS5-brane charge. 

The holomorphic functions ( |5.4j ) are supposed to be considered as transition functions 
modifying the complex contact structure of the twistor space of the c-map. 5 Our results 
suggest that, up to a universal factor t P//?l , they can be identified with (the quasiclassical limit 
of) the B acker- Akhiezer functions, provided the Planck constant and the NS5-brane charge 
are related as 

fT l = 8nk. (5.7) 

Amazingly, a similar relation between a quantization parameter and the NS5-brane charge 
has been suggested in |L9| on the basis of a beautiful connection between the geometric 



quantization of cluster varieties f55[l , motivic wall-crossing formula j|0], ^TJ, and complex 
contact structure. 

Let us recall also that in MQM the two Backer-Akhiezer functions can be seen as per- 
turbed one-fermion wave functions in the two chiral representations where x + and x_ appear 
as diagonal operators, respectively. Thus, the wave function in the right representation en- 
codes the effects due to NS5-branes of charge k and the (conjugate) wave function in the left 
representation does the same for anti-branes. The precise relations read as follows 

where the complex conjugation does not act on the fields and at the end one should take 
^mqm = (87ri/c) _1 (see footnote [I]). This identification is in the beautiful agreement with a 
general expectation that the NS5-brane partition function should behave as a wave function 
42, 43, 44, 45] and, furthermore, suggests that NS5-branes might have a fermionic nature. 

Moreover, the duality with MQM provides an interesting interpretation for the contours 
C±. In [^J they have been chosen by hand to produce the correct instanton effects ( |5.6| ) 
and their meaning was completely unclear. On the other hand, in terms of the MQM chiral 
coordinates, these contours turn out to be very simple: they join x ± = oo to x^ = 0. The 
starting points correspond to the asymptotics of the Fermi sea, whereas the end points are 
nothing else but the intersection points of the Fermi level with the lines of the perturbative 
stability. Namely, the inverse oscillator potential V — — \x 2 has two infinite "wells". Ignoring 
non-perturbative effects, the movement of a particle is confined to one of them provided its 
initial conditions satisfy x 2 — p 2 = 2x + x_ > 0. Thus, the contours C± correspond to the parts 
of the Fermi profile which bounds the fermions leaking to the other side of the potential. We 
demonstrated this situation on Fig. |[ 

new feature here comparing to what we presented above is that these transition functions are not 
associated to an overlap of any two open patches. In that case they would be integrated over a closed contour 
surrounding one of the two open patches. Instead, the contours C± are open. Such contours typically arise 
when the Darboux coordinates develop branch cut singularities. In particular, this is the case for the twistor 
space description of D-brane instanton corrections to the local c-map |5(| and BPS instanton corrections to 
the rigid c-map (s). 
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Figure 2: The movement of free fermions in the inverse oscillator potential. The second picture 
shows the profile of the Fermi sea on the two sides of the potential for the choice of the couplings 
t±i = =Fl. t± are the points of intersection with the lines = 0. They bound the parts of the 
Fermi sea (marked by a darker color), with the phase space coordinates satisfying x + x_ < 0, which 
leave to (come from) the opposite side of the potential. This is illustrated on the third picture which 
shows the time evolution of the Fermi sea from the second picture. One can see how a part of the 
sea leaks from the right to the left. 



5.1 Inclusion of the one-loop correction 



It is very easy to include in the above description the one-loop ^-correction to the universal 



hypermultiplet geometry [|7|]. As was shown in [28|, in the twistor formalism it appears 
as a non-vanishing anomalous dimension c, which is fixed by the Euler characteristic of the 
Calabi-Yau, c = — This anomalous dimension has only two effects: it adds a logarithmic 
term to a (see ( ft.Tp ) and shifts the contact potential by a constant. As discussed below 



( |3.8|) , this amounts to shifting in all above expressions p by the constant c. In particular, the 
one-loop corrected metric corresponds to the solution of the Toda equation T = log(p + c) 
and therefore leads to a modification of our dictionary by setting p, — p + c. The only place 
where p remains unshifted is the relation ( |5.3| ) between a and the WKB phases since the shift 
is canceled by the logarithmic correction to a. 

The relation between the NS5-brane instantons and the Baker- Akhiezer function works 
in a similar way as above. One should only take into account that the function ( |5.4|) should 
be multiplied by (£ =F £) 87Tck [27]. As a result, the relation ( |5.8|) generalizes to 



H, 



NS5 
k,+ 



-nika 



(5.9) 



6. Discussion 

In this note we proposed a correspondence between the QK space given by the local c-map 
and the matrix model description of the c = 1 string theory with the sine-Liouville pertur- 
bation compactified at the self-dual radius. This duality is based on a dictionary between 
the twistorial description of four- dimensional quaternionic spaces having an isometry and the 
Lax formulation of the dispersionless Toda hierarchy. Under this duality, the coordinates 
parametrizing the quaternionic space appear as couplings of the c = 1 string theory, whereas 
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the metric is encoded in the r-function of the hierarchy. Besides, the coordinate parametriz- 
ing the CP 1 fiber of the twistor space appears on the matrix model side as a uniformization 
parameter of the Riemann surface given by a complexification of the profile of the Fermi sea 
build from the free fermions of the matrix model. 

While this curious duality is very interesting, it would become useful if it can be applied 
to extract some information about physical effects on one side from the known results on the 
other side. In particular, it is very important to understand whether it can be extended to 
include quantum effects. While the inclusion of the one-loop correction to the c-map is almost 
trivial, we also provided an identification of the one-instanton corrections due to NS5-branes 
with the Baker- Akhiezer function of the Toda hierarchy, which is realized in the c = 1 string 
theory as the wave function of a free fermion in the chiral representation. This relation points 
towards a fermionic interpretation of NS5-branes which might have important implications 
for the study of their dynamics. In particular, one might expect that multi-instanton effects 
should be expressible in terms of a Slater determinant. 

On the other hand, it is not clear how D-brane instantons to the hypermultiplet geometry 
can fit into this correspondence (see [[32|] for a review of the current understanding of the non- 
perturbative effects in the hypermultiplet sector). The problem is that there are two string 
couplings in the game which have different scaling. One is the coupling constant of the c = 1 
string theory which scales as g^ =1 ~ ^ l . The second string coupling describes quantum 
corrections to the hypermultiplet moduli space and is set by the dilaton to g s ~ p -1 / 2 . 
Recalling that fi = p, one obtains that, if D-instanton corrections to the c-map have a 
counterpart, it should scale as e , which seems to be rather strange. 

Moreover, one of the key ingredients determining the structure of instanton corrections 
on the hypermultiplet moduli space is the Heisenberg symmetry. In particular, it requires 
that the coordinate z, describing the RR-fields, lives on a torus. Our correspondence then 
implies that the sine-Liouville couplings t±\ should be periodic. It is not clear at all how such 
periodicity condition can arise in the c = 1 string theory. 

Thus, if the duality can be extended to include quantum effects, it should involve a non- 
trivial deformation of the c = 1 string and may even require to replace it by a different 
physical system. Some indications towards this conclusion can be seen already here. First of 
all, note that the NS5-brane instantons ( |5.6|) have the same strength as the tunneling effects 
in the free fermionic system given by MQM. Due to this, one may expect that the fermions 
should fill both "wells" of the potential, which leads to consideration of the supersymmetric 



OB two-dimensional string theory |j48|| . Second, there is a problem that the c = 1 string and 
the c-map impose different reality conditions on the Toda times: t* ±l = t±\ and t\ = 
respectively. This hints that the proper system dual to the c-map is not the c = 1 string, but 
the so-called Normal Matrix Model. It is given by the same solution of the Toda hierarchy, 
but describes a different real section of the same holomorphic structure [^, [50| . 

It is important also to understand whether the presented results can be extended beyond 
the restriction to four dimensions. The main complication appearing in higher dimensions is 
that quaternionic geometries do not coincide anymore with self-dual Einstein spaces. Never- 
theless, one may hope that the Toda hierarchy still describes them upon switching on higher 
times t±k- Even if this is not the case, there is still a possibility that an extension of this 
hierarchy to a more general one will do the job. 
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So far we were looking for an interpretation of the physical effects in the hypermultiplet 
sector of type II string theory from the point of view of the c = 1 string or its MQM realization. 
One can consider also the opposite problem. In particular, in the c = 1 string theory with the 
sine-Liouville perturbation the non-perturbative corrections have been studied in |51, [52|. 6 



One of the nice aspects of the resulting picture is that the instantons are associated with 
singularities of the Riemann surface already mentioned above [54. 55]. It would be interesting 



to see whether such geometric interpretation can be translated to the twistor space of the 
non-perturbative moduli spaces in N = 2 theories. 

Another issue to be understood is that on the Toda or c = 1 string side, the presented 
duality involves so far only the dispersionless limit. As was reviewed in section |2.1| , there is 
a natural quantization of this classical structure generating a perturbative expansion in h. 
In c = 1 string theory the H- correct ions are interpreted as contributions from worldsheets 
of higher genera. A natural question is: Is there a meaning of a similar deformation for 
quaternionic geometries? 

There are several indications that this deformation corresponds to the inclusion of NS5- 
brane instantons. Indeed, the Heisenberg symmetry mentioned above, in the presence of 



a non-vanishing NS5-brane charge, becomes equivalent to the quantum torus algebra [ [45 
Therefore, one can expect that switching on this charge corresponds to promoting certain 
structures from classical to quantum and this quantization is equivalent to the one in the 
Toda hierarchy. This expectation is strengthened by the observation that our duality implies 
the relation ( |5.7| ) between the NS5-brane charge and the Planck constant. Note also that the 
Toda hierarchy has been shown to incorporate the quantum torus algebra in several physically 



relevant situations 
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